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1 Introduction 



In the previous paper referred here as I, the contribution of the core nucleons to the anapole 
moment (AM), arising from the core polarization effects, has been calculated in the random- 
phase approximation (RPA). Recently, the first measurement of the AM of Cesium has been 
reported 0. The immediate application of this measurement was the attempt to deduce the 
pion-nucleon weak parity- non-conserving (PNC) coupling constant [^], ||. The comparison 
of the measured anapole moment value with the one calculated using a pure single-particle 
model leads to a value for f n that exceeds by a factor 4 the value deduced from a parity 
violating measurement in 18 F |J. More sophisticated comparison was made in Ref. |J, where 
the results from 18 F, the AM of 133 Cs, and the upper bound for the AM of 205 Tl J7J were used 
to deduce both and f p PNC coupling constants. The combination of the coupling constants 
was found which satisfies both 18 F and 133 Cs experiments and which is barely in agreement 
with theory. These values, however, are inconsistent with the constraint obtained from 205 TZ 
measurement. This situation, even independently on 18 F experiments, raises the question how 
accurate is the theory of nuclear anapole moments. Here we address this question and present 
more accurate calculation of the many-body effects. 

There are different contributions to the value of nuclear AM. The most significant and the 
least model dependent contribution comes from the valence nucleon. It is stable under variations 
of the single-particle potential ||. Within the single-particle model it is often convenient to 
use the leading approximation (LA) || producing the simple analytical expression for the 
correction to the single-particle wave function Sip due to the PNC potential. The LA was used 
in I in calculations of the core nucleons contribution to the anapole moment. This is another 
contribution which is not negligible. The accuracy of the LA was estimated to be within 20% 
in ||. We found, that although this is true for the upper level of the spin-orbit doublet, for 
the lower level the difference between LA and the exact correction 5ip can be significant. This 
demands improving the calculation of the core nucleons contribution. 

Another phenomenon that has to be accounted for is the pairing. The reason for this is that 
for 133 Cs nucleus there are single-particle states very close to the Fermi surface. The transitions 
to these states produce large contribution to the core polarization. In the presence of pairing 
these transitions will be strongly reduced, thus affecting the value of the polarization effects 
and the nuclear anapole moment. 

Finally, we discuss the constants of the short range effective PNC interaction used in nuclear 
structure calculations. These constants were obtained from the finite range meson exchange 
forces by comparing typical matrix elements both finite range and short range interactions 



1C| |11] , |12fl - The factors W p ~ 0.4 and W n ~ 0.16 were found for tt- and p-exchange forces. 
These values, however, were obtained for a-particle and the extension to heavier nuclei should 
be checked separately. We found using the same procedure that the factor W n is state dependent 
and it is different for 133 Cs and 205 Tl. 

The paper is organized as follows. First, we discuss the accuracy of the LA. Then, we present 
the exact calculations of the core polarization contribution and demonstrate the necessity to 
include the pairing. Finally, we discuss the constants of the short range PNC interaction. 
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2 Core polarization contribution 



For completeness of the discussion let us to remind the basic equations from I. There are 
several contributions to the AM arising from different parts of the electromagnetic current. 
Apart from magnetization current which gives the main contribution, there are contributions 
from the convection current, the spin-orbit current and the contact current arising from the 
momentum dependence of the corresponding parts of the Hamiltonian M : 
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Here { , } is anticommutator, < > is the proton mean squared radius, p is the central 
nuclear density, f(r) = p(r)/p - nuclear density profile, and Uf s n = Ujf = 134 Me V ■ fm 5 is 
the proton-neutron constant of the effective spin-orbit residual interaction ||. Gf is the weak 
interaction Fermi coupling constant, m is the proton mass, and fi a is the nucleon magnetic 
moment. The contact current contribution arises from the velocity dependence of the effective 



nucleon-nucleon PNC forces || 14 1 



The operators in Eq. (fj) have 3 types of angular dependence. All of them, except the 
convection ones, are spin dependent. All of them except a c are of El type. The operators 
created by the contact current are of Ml type. The convection terms in Eq. (0) differ from 
the one used in I and ||. Here we accounted for the recoil correction by subtracting the part 
related to the motion of a nucleus whole. 

As in I, we use for the radial dependence of the effective operators related to Eq. ([!]) 
the notation V[aj](r). The RPA equations for the effective radial fields V[aj](r) are (see Eqs. 
(17), (18) in I) 

V = V ([ai])+F a AV, (2) 
8V = F W AV + F s 5AV + F s A5V. (3) 
Here F w is the PNC nucleon-nucleon interaction 

Q I 

Fw = "71 4m ^ ^ 9ab(Ta ~ 9baCTb > ' ( Pa ~~ P6 )' S ( Ya ~ r& )} 

+ <4K X <r b ] ■ V5(r ffi -r 6 )), (4) 
and F s is the residual spin-spin interaction 

F s (ab) = C (g + g' T a • r b ) er a ■ cr b 5(r a - r b ). (5) 
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As in I, we use C = 300 MeV ■ fm 3 , go = 0.6, and g' Q = 1. Averaging the interaction Eq.(^) 
over core particles we obtain a single-particle PNC potential 

^W = ^«-P)./M}, (6) 

where po is the central nuclear density, f(r) is the nuclear density profile and the constants g a 
are related to the interaction constants g a b (see Eq. (4) in I). In the absence of pairing, the 
static particle-hole propagator A(r, r') has the form 

A(r,r') = E^(r)^(r')^^^(r)^(r), (7) 

vv' tv €v ' 

where n u are the occupation numbers, e u and ?/v(r) are the energies and the wave functions of 
the single-particle levels. The Feynman diagrams corresponding to Eq.(fJ,[|) are shown in Fig.l. 
Eq.([J), which is shown in the top diagrams (I) in Fig.l, describes the usual RPA renormalization 




(a) (b) (b') (c) 

Figure 1: Feynman diagrams corresponding to Eqs.(2),(3). 

I. The anapole moment renormalization diagram describing the core polarization effects. The 
filled circle is the bare anapole operator. The open triangle is the dressed one. The open square 
is the spin-spin residual strong interaction. 

II. Additional contribution to the anapole moment due to parity violation in the core states, a) 
- the direct contribution of the two-particle effective PNC interaction, b) and b') - the effective 
contribution due to parity violation in the core states, c) - the renormalization diagram. The 
shaded circle is the single-particle PNC interaction. The shaded square is the two-particle 
effective PNC interaction. 

of the bare operator Vo[aj](r). The next Eq.(|3]), which is shown in the bottom diagrams (II) in 
Fig.l, describes an additional contribution from the core nucleons arising both from the direct 
P-odd nucleon-nucleon interaction F w , see Fig.l (II a), and from the P-even residual interaction 
via admixture of opposite parity states to the wave functions of the core nucleons,see Fig.l (II 
b), (II b'). The last term in Eq.(|3|), see Fig.l (II c), is responsible for the renormalization of 
these contributions. 
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It is worth mentioning that the correction 5V has the parity opposite to V. The change 
in parity happens because SV is created from V by the interaction that does not conserve 
parity While AM is El type operator, the correction 5V is Ml type. For this reason the 
renormalization due to the core polarization will be different for V and SV, since different 
transitions are involved in the kernels of integral equations (Q) and @. Eq.(H) describes also 
the renormalization of the operators a c . They don't have SV since they themselves are of first 
order in the weak PNC interaction. The AM value is given by the sum of all these terms: 

a, = (5^\V[ai]\^) + (VIVMISV) + (V|W[ai]|V>, ( 8 ) 
where the index i means here s, conv or Is. For the contact contribution we have 

a c = (V>|V[a c ]|V>>. 

2.1 Leading approximation accuracy 

Eq.flU) was solved in I using LA in calculations of 8 A. The accuracy of LA was estimated in || 
within 20%. This estimate was maid for odd valence nucleons in a set of nuclei. Occasionally, 
all the valence nucleon levels under discussion were the upper levels of the spin-orbit doublet. 
The typical difference between exact SR(r) and rR(r) in LA for the upper level of the spin-orbit 
doublet is shown in the left plot of Fig.2. However, the sum in 5 A goes over all states. For the 




2 4 6 8 10 12 14 16 ' 2 4 6 8 10 12 14 16 



I (fm) r (f m ) 

Figure 2: Admixture of the opposite parity 5R(r) for different components of the spin-orbit 
doublets. Left panel - 2d3/2 level, right panel - lgg/2 level. LA is shown by the dashed lines. 
Full lines show the exact solution. 

lower states of the doublet the difference in the peak heights can reach a factor 2 as one can 
see in the right plot in Fig.2. In all cases the exact correction and the LA are peaked near the 
surface. However, the ratio of the peak heights differs considerably for the upper and the lower 
states of spin-orbit doublets. 

The systematic study of this ratio is presented in Fig. 3. Here we plotted the ratio of the 
peak heights of the exact SR(r) and rR(r) in the LA for a set of proton states in the Woods- 
Saxon potential. While for the states with j = 1 — 1/2 the ratio remains close to 1, for the states 
with j = 1 + 1/2 the ratio goes down with increasing orbital angular momentum I. For the 
state l/in/2 the exact correction 5R(r) differs more than by factor 2 from the correction rR(r) 
calculated using LA. From these results we see that the estimated accuracy 20% in calculation 
of S A cited in I is too optimistic and more accurate calculation is necessary. 
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Figure 3: Ratio £ of the maximum in exact 5R(r) to the maximum in 5R(r) found in LA for a 
set of the proton single-particle orbitals. 



2.2 Exact 5 A 

Let us start with the expression F S AV in Eq.(|2]). In coordinate space it can be presented as 
follows 

{F s AV[m\) b {r) = Cg£<TX c (r), 
where C and g bc are the parameters of the residual interaction (13), and 



(9) 



In the following derivation we keep the external frequency uo. It will be put to zero in the final 
result. Introducing the single-particle Green's function 



G(r r'\E) = £ 



E-e„ 



we can rewrite the Eq.(^) in the following form 

X(r) = 5>„ [ d 3 r {^(r)ffG(r r'le, + u)V{v')^ v {v') 

+ xPt(r')V(r')G(r' r\e u - w)tnp v (r)} . (10) 

The correction to Eq.(pH|) arising from the single-particle PNC potential can be obtained 
using the correction to the single-particle wave function ^v(r) defined in I: 



where 



5i/) v (r) = -t£ a (ern)5R u (r)Q u (n) = i£ a 5 R u (r)Qv(n) 



e G 

<ia — —^9aP0- 



(11) 
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We used here the relation 

% m ( n ) = -(cm)^ 7m (n) 
where I = 2j — I. For the correction <5X(r) we find then 



5X(r) = £Xi(r) + £K 2 (r) 



where 



and 



vv' 

(fit(n)fffi I y(n))G? I /(r,r / |e I/ + a;) < z/|V(r')|z> > 5R v (r')R v (r) 

- {nl{n)<Tn v ,(n))G„,{r,r'\e v + u) < v'\V{r')\v > 5R u {r)R u {r') 
+ 5R v (r)R v (r') < u\V(r')\u' > (nl(n)an p (n))G v ,(r'r\e v -w) 

- R„{r)5R„{r') < v\V(r')\v' > (fit / (n)«7fi 1/ (n))Gf I/ /(r / r|e v - w)} , (12) 

5X 2 (r) / rfVj^W^^r'l^ + ^y^)^^) 

+ ^(rWO^* r'K - w)^(r)} . (13) 



The angular brackets in Eq.fll2|) denote the matrix elements in angular and spin variables 
only. In the following calculations for spherical nuclei, it is convenient to separate the angular 
variables by expanding both 5Xi(r) and 5X 2 (r) in a set of the vector spherical harmonics 
Y&,(n), 

5X i (v) = J2xt(r)Y5 M (n). (14) 

L 

In Eq. (fn|) the total angular momentum J is fixed, J = 1, and the sum goes over L only. 
Multiplying Eqs. (H -HI) by YjJ^n), integrating over angles and separating the dependence 
on all angular momentum projections M via Wigner-Eckart theorem, we obtain for the reduced 
functions of r 

vv 1 

x (G^(r,r'|e^ + ^) + G^(r,r'|e^-cj)). (15) 
Here, double vertical lines mean the reduced angular matrix elements. The tensor operators 

Tj M = {a®Y L } JM (16) 



Tj M are defined as in I: 



Similar expression can be obtained for X 2 (r) 



x / r'W(5G^,(ry|e, + o;) + ^^ (17) 
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In Eqs. ( |T5| , |TTD 5G is the correction to the Green's function of the Schrodinger equation for 
single-particle orbitals due to the PNC potential Eq.(EI). In first order the correction is 



5G(r,r') = J d 3 iG(r,x)lf(x)G(x,r') 



Separating the angular dependence we obtain for the radial correction 

( <— > \ 
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Here k = (I — j)(2j + 1). For obvious reason the, PNC potential couples only the states with 
the angular momenta I and I. 

The final equations for the radial functions v [a*](r) and dvL[a-i](r) are as follows 



v°[a<](r) - J2 9? / T^drM^rVkKO = «(r); 
. — JO 

b=p,n 

KNW-EJo'E / r' 2 dr'^ L ,(r,r')K'N(0=K a W(r), 

h—„ „ r i—n o J 



(19) 
(20) 



b=p,n L'=0,2 

where the radial particle- hole propagator A b LL , (r, r') is 



A b L Jr,r' 



\C E n) ln m\T^\\fl')m\T^ \\j'irRMr)R jln (r') 



O " Z — ' J 

° jlnj'V 



x (Gfi'far'-^jin + u) + G f i> (r,r'; e JJn -w)), (21) 

with u; — > 0. Here C is the normalization constant of the residual interaction Eq.(^) and 
Gji(r,r'; e) is Green's function of the radial Schrodinger equation. 

Eqs.(r^|),(PD|) are similar for all anapole operators Eq. ([]]), except for the convection current 
operator which will be treated separately. The right-hand-side of Eq. ( |I9|) for the spin current 
anapole operator is 

<[a,](r) = --^^- 



-r. 



rn 



(22) 



The transition from the Cartesian vector product in Eq. (Jl|) to the tensor operator T± M results 
in the factor —iyj8n/3. The factor is common for the spin and spin-orbit anapole operators. 
The radial dependence follows Eq. (|]) . The right-hand-side of Eq. (|20"D is given by 



b=p,n 



(23) 



In Eq.(P3"|) we omit the direct contribution of the two-particle effective PNC interaction 
corresponding to the diagram (Ha) in Fig.l. According to the estimates made in Ref.fT3|| this 
contribution is small, it does not have the factor A 2 / 3 . The estimate in [13j was obtained with 
harmonic oscillator wave functions. We calculated the diagram (Ha) for y°[a s ] from Eq.(|l]) 
with the Woods-Saxon wave functions. For 133 Cs nucleus we obtained the contribution 



8k Fw 



AAg pp + 0Ag pn + 4.9g np - 1.6g' 
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-0.015 



for the "best values" of the coupling constants |I3. This number should be compared with 



the single-particle value of k s = 0.32, and it gives the contribution less than 5%. Similar 
calculations give 5k Fw = 0.003 for 205 TZ, 5k Fw = 0.023 for 209 Bi, and 8k Fw = 0.0085 for 207 Pb. 
All these values are in agreement with the estimates of Ref. [|1^ . 

The convection current anapole operator Eq.(|l|) expressed via tensor operators has the form 



a mnv)M 




-^ +1 r y -i <r > }YlM )- (24) 
= < rl > (1 { .,y !m } - ^ ( | + . 

The operators Eq. (p4|) are non-local and do not include spin. For the last reason they are 
not renormalized directly by the residual interaction Eq.(^). However, due to the spin-orbit 
potential a new spin dependent contribution to AM is generated in first order in the residual 
interaction Eq.(^). Its tensor structure is given by the same tensor operator as that for the spin 
current AM a s . The radial dependence of this polarization contribution is given by 



b=p,n 



/■oo 

x / r'^rX^O^IIaLJkO^KH^ + ^ + ^^r'le,-^)), (25) 
Jo 

where uj — > 0. This operator undergoes the usual renormalization described by Eq.(plJ) pro- 
ducing in the next orders in the residual interaction the renormalized radial operator Vp 0l (r). 
Although the direct contribution of the convection current to AM is small, the polarization con- 
tribution Eq.(^) is larger and has to be included. Total operator generated by the convection 
current has the following structure 

v b Won V ]{r) = a b conv + v b pol {r)Tl (26) 



where a conv is given by Eq.(Bl) and v b ol (r) is the dressed polarization operator Eq. (551) . The 



tensor operator is defined by Eq.(|l^) and we omit the total momentum projection. This 
very operator Eq.fl^Tf) should be used in the Eq.(P3"|) for the right-hand-side of Eq. (|2H|) . 

The values of the AM resulting from the solutions of Eqs.([19D, (p0|) are summarized in Table 
1 and Table 2 for 133 Cs and 205 TZ nuclei. For comparison, in the row labeled by I we listed the 
results of previous calculation with the use of LA. The results of present calculations are listed 
in the row labeled by II. 

The value cited in column V — Vq is the renormalization of the single-particle (column s.p.) 
value due to the core polarization by Eq.(|3|). One can see that the polarization contribution is 
about half of the single-particle one and it has an opposite sign compared to the single-particle 
contribution for all spin dependent operators. This is in accordance with the repulsive nature 
of the spin-spin residual interaction Eq.(|]). In this case the core produces a screening of the 
valence nucleon spin. 

The contribution of the convection current in the Tables 1 and 2 was listed together with 
the polarization effects discussed above. The single-particle value is small compared to the 
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magnetization current contribution. Let us note that it differs from the corresponding value 
cited in I. The difference, although small, comes from the center of-mass-motion that was 
not excluded in I. It is interesting to note that for the convection current, the polarization 
contribution has the same sign as the single-particle one. This is in contrast to the spin case 
where the polarization contribution has always the opposite sign compared to the single-particle 
value. The difference between these two cases lies in the polarization loop which is non-diagonal 
for the convection current case. It has the spin operator in one vertex and the convection part 
of the AM (see Eq.(|IJ)) in the other. For such non-diagonal loop the sign is not fixed. 

The main conclusion that can be drawn from the present results is that the LA overestimates 
the contribution of 5V approximately by a factor 2. The summed value K to t has decreased by 
50%, as compared to its single-particle value. Thus, the effect of the core polarization is 
considerable. 



3 Pairing effects 

In general, the pairing effects are important only for the transitions near the Fermi surface. 
Therefore, they can be neglected in calculations of 5ip and in the polarization loop in Eq.(|T9|) 
where the transitions to at least next shell are involved. The situation is different for Eq . (|20|) 
where the transitions within the open shell are allowed. Good example is 133 Cs where the 
polarization loop includes the transitions between the states with the excitation energy about 
one hundred KeV. Transitions to such close levels contribute significantly into polarization 
loop creating large response in 5V. These transitions will apparently be suppressed by pairing 
correlations. For these reason the polarization loop in Eq.(pO|) should be modified in order to 
include pairing. 

The modification is straightforward and can be done similarly to jn|. The particle- hole 
propagator for a T-odd channel in the presence of pairing has the form 

A((r T f ) = -y2 ^ UvVv ' ~ ^^) 2 ^( r )^ , ( r )^ , ( r/ )^( r/ ) (27) 

Ev + E v > 



where u v and v u are the Bogolyubov factors, and E v = J (e v — ji) 2 + A 2 where A and \x are 
the pairing gap and the chemical potential. Below, we shall neglect slight state dependence of 
the pairing gap and put A = const. Let us divide the single-particle space into 3 regions. Let 
P v be a projection operator onto the region near the Fermi surface, where pairing produces 
a significant effect. It is convenient to include here all bound states above the Fermi surface. 
Let Q v be a projection operator onto the hole states below the pairing region. In this region 
we neglect the pairing effects. The states in continuum will be projected by 1 — P v — Q v . 
Introducing the identity 

1 = (P V + Q V + 1-P„- Qv)(Pv> + Q V > + 1- P v > - Qy>) 
into Eq. (|27|) we obtain the expression for the particle-hole loop as a sum of 7 following terms 

A(r,r') = -E^(r)^(r)^(rO^(rO x [ ^'^f 



ul _ ul, 



PvQv' + r= ; P u 'Qv 



E u - e u > + fi E v > - e u + fi 
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vl „ ,„ „ „ s . vl 



+^r~^ Pu{l - Pv< - Qu>) + ^—^ Ps(l -P v - Qu) 

E v + e v > - n E v > + e v - ji 

1 -Qsil-Py-Qu) — QuiX-Pw-Qv^Y (21 



In Eq.([28|) we put A = in all terms containing Q v . The sum over the whole single-particle 
space is present in the last 4 terms containing the unit projection operator. Combining these 
terms together we can present Eq.([28|) as a sum of three different terms 

A(r, r') = Aq(t, r') + A P (r, r') + A e (r, r'), 

where 

Aq(t, r') = £ QAM^Kr'Mr', r|e„) + G(r, r'|e,)^(r')^(r)} (29) 

V 

includes the transitions from the deep holes, 

A P (r, !•') = -£ P^{^(r)^(r')G(r', r|/i - £„) + G(r, r'|// - ^)^(r')^(r)} (30) 

includes the transitions from the pairing region, and 

A e {r, r') = -E ^(r)^(r')6^^(r')^(r) (31) 

with 

+ b ^ - TT-^ + 

r vl ' 1 1 ]P„Q„. 



— £v + + e ^ ~ A 1 e !/ — e */ 

includes the transitions from the pairing regions to the region of deep holes. In all these terms 
the sum goes over a finite range of the single-particle space and can be calculated directly. 

Strictly speaking, in the presence of pairing the particle-particle channel should be added 
to Eq. ([Z5D and the equation for V has to be modified. With the particle-particle channel a 
new interaction appears and the total number of equations doubles. The parameters of 
the interaction F^ can be found from masses of the nuclei that differ by two protons or two 
neutrons |T3. This procedure, however, determines the spin-independent part of the interaction 



F^. The spin- independent interaction in the short-range approximation does not renormalize 
the spin-dependent operators. As for the spin-dependent part of F*, its magnitude remains 
unknown up to now and for this reason we did not include it in our equations. 

For the right-hand-side of Eq.fl2"3"|) one can obtain 5 A from Eqs.(P9"D-(|5T]) using in first order 
5ip by Eq.fllT]) and 5G by Eq.(|l|). 

In Fig. 4 we plot the difference v%(r) — Vq s (t) for the spin part of the AM operator both in 
the absence and in the presence of pairing. As expected, the difference between these two cases 
is insignificant. For SV the situation is somewhat different. In Fig. 5 we plot the right-hand-side 
of Eq.(|23|) for the transitions with L = (left plot). In this case the difference is also small. 
However, for the transitions with L = 2 (right plot) the difference is significant. This is the 
direct influence of pairing that reduces the transitions to closely lying levels. The final value of 
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Figure 4: The difference v p s (r) — Vq s (t) for the spin part of the AM operator in the presence of 
pairing (solid line), and in the absence of pairing (dashed line). 




Figure 5: The contributions to proton X L [t>o s ] with different L in the absence (dashed line), and 
the presence (solid line) of pairing. The interaction constant g p put equal to 1 for illustrative 
purpose. On the left plot L = 0, on the right plot L = 2. 
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AM does not changed strongly because the contribution of X L [t>g s ] with L = 2 to AM is small. 



The results of the anapole moment calculations with pairing are listed in the Table 1 for 
133 Cs nucleus in the row labeled III. For the pairing gap A we used the values A p = 0.95MeV 



for protons and A n = l.ObMeV for neutrons from Ref.||17||. As expected, the pairing does not 
influence much the matrix element of the anapole moment operator V. For SV the effects are 
relatively larger. As mentioned above, the pairing reduces the transitions near the Fermi surface 
with small AE. For the Cs nucleus this happens only for the proton transition 1(77/2 — 2g? 5 / 2 , 
Al = 2 that does not contribute significantly to the anapole moment. 



4 Parameters of the PNC nuclear forces 

The constants g a b and g' ab of the PNC interaction (f|) should be, strictly speaking, treated as 
phenomenological ones and should be found from experimental data. We can, however, try to 



relate them to the parameters of the free nucleon-nucleon interaction |fL8[ . This relationship, 
as obtained in is: 

g pp = -(ji + 2)W p A p ti± 



9pp — 9nn — 9nn ~ 9pp 



g pn = -(2fi + l)W p A p h° p + W^U (32) 
g np = -(2/i + l)W p A p h° p - W W A W U 
9' P n = 9n P = (V - VWpAph , 



where 

A A 9 * 



p G F m 2 p G F m\ 

are the dimensional constants, /i is the isovector nucleon magnetic moment, and h° p , and f w are 
the PNC rho-nucleon and pion-nucleon couplings. The dimensionless factors W p and were 
introduced to normalize the matrix elements of the interaction Eq.(§) to the matrix elements 
of the original finite range interaction ||18|| . 

Following Ref.||, in Fig. 6 we plotted the extracted values of coupling constants H\ and 
Hp in the notations of Adelberger and Haxton ||. They are related to /„- and h Q by 

rrl _ 9 tt fit tjO 9p^p 

n ~ a/32' p ~ 2 ' 



where g n and g p are the strong coupling constants. The bands corresponding to the 205 TZ 
and 133 Cs data are slightly different from those extracted in Ref. ||. There, our previous 
calculations M were used to extract the coupling constants from the anapole moment data. 

Our improved calculations did not change the general situation. The coupling constants 
extracted from 133 Cs and 205 77 data still look inconsistent independently on 19 F data. This 
situation eventually raises a question how reliable is the theory of nuclear anapole moment. 
Our calculation has been done using the random-phase approximation. There are two more 



calculations accounting for many-body effects [20|, [21] within the shell model approach. In 
Ref. [|2l| the shell model basis used for calculation of the anapole moment of 205 TZ was large 



enough to account simultaneously both for the single-particle AM and the core polarization 
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Figure 6: Weak coupling constants H\ and H® extracted from 18 F (light band), 205 Tl (medium 
band), and 133 Cs (dark band) experiments. Also shown are the DDH JT8J "best" values (square) 
and "reasonable range" (box) 



effects. The value obtained in Ref. [21 for the spin part of the anapole moment of 205 Tl is 
k s = 0.35. Our calculation, where we use completely different quasi-particle interaction and 
the RPA, gives k s = 0.37 for the "best" values of the coupling constants. Such close values 
obtained in completely different approaches give hope for a weak model dependence of the AM, 
although we cannot exclude that this is just the coincidence. 

The relation between effective PNC interaction constants and meson-nucleon coupling con- 
stants given by Eq.(^) is less reliable. The normalization factors W n and W p were determined 
by comparing matrix elements of the interaction Eq.(^j) and the finite-range DDH interaction 
of Ref. |L8| for N — 4 He scattering at low energy. Therefore, we hardly can expect them to be 



constant in a broad range of all bound single-particle states. In order to check the accuracy of 
this procedure we calculated as a ratio of typical matrix elements of the finite range DDH 



181 interaction and the zero range interaction Eq.(£|). The results are shown in Fig. 7 where the 
normalization factor W n is plotted as a function of a single-particle energy for the proton and 
neutron states. The factor W„ is really state dependent. Only near Fermi surface W n = 0.16. 
For the states with lower energy it becomes larger reaching the value ~ 0.4 — 0.5 for lsi/2 
states. For this reason one should not use such a simple relation as Eq. ( |32f) for determination of 
the PNC meson-nucleon coupling constants. There is an additional reason why Eq. (pSf) should 
not be used to extract the interaction constants. As it was shown in p2[ ^3[ strong renor- 
malization of the PNC nucleon-nucleon interaction can exist in nuclear media. As a result, 
the neutron PNC potential constant g n (which is small when estimated using Eq.(|32"D) can be 



comparable to the proton constant g p . (See also discussion of the subject in Ref.f24fl.) From this 



point of view, the measurements sensitive to g n would be extremely interesting. These could 
be the measurements of the anapole moment of a nucleus with an odd neutron or, another 
possibility, the measurement of the neutron spin rotation in helium |TD . 



In order to obtain the measured value of the AM of 133 Cs the PNC interaction constants 
Eq.(|^) should be increased approximately by a factor of 2 as compared to their Eq.(|3"2"D "best 
values". It is interesting to note that a similar conclusion has been obtained from statistical 
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Figure 7: Normalization factor W n as a function of bound state energy. Proton states are 
connected by dashes. Neutron states are connected by dots. 



analysis of the PNC effects in compound nuclei |25 . 

In Table 3 we show the summarized results for a set of the proton odd and the neutron odd 
nuclei. Here we list just the sum of all contributions. For Ba isotopes the calculations were 
performed including the pairing. It is worth noting that for nuclei with an odd neutron due to 
the core polarization some contribution proportional to the proton coupling g p appears in the 
anapole moment. In case of small g n it can be of the same order as the direct g n contribution. 

Acknowledgments 

The discussions with LB. Khriplovich are greatly appreciated. One of the authors (VD) ac- 
knowledges the discussions with B.A. Brown and V. G. Zelevinsky. This work was supported 
by the RFBR grant 98-02-17797. 

References 

[1] V.F. Dmitriev, V.B. Telitsin, Nucl. Phys.A 613 (1997) 237. 
[2] C.S. Wood et al, Science 275 (1997) 1759. 
[3] W.C. Haxton, Science 275 (1997) 1753. 

[4] V.V. Flambaum and D.W. Murray, Phys. Rev. C56 (1997) 1641. 

[5] E.G. Adelberger and W.C. Haxton, Ann. Rev. Nucl. Part. Sci. 35 (1985) 501. 

[6] W.S. Wilburn and J.D. Bowman, Phys. Rev. C57 (1998) 3425. 



15 



[7] P.A. Vetter, D.M. Meekhof, P.K. Majumder, S.K. Lamoreaux, and E.N. Fortson, Phys. 
Rev. Lett. 74 (1995) 2658. 

[8] V.F. Dmitriev, LB. Khriplovich and V.B. Telitsin, Nucl. Phys. A577 (1994) 691. 

[9] V.V. Flambaum and LB. Khriplovich, Zh. Eksp. Teor. Fiz. 79 (1980) 1656. 

[10] V.F. Dmitriev, V.V. Flambaum, O.P Sushkov, and V.B. Telitsin, Phys. Lett. B125 (1983) 
1. 

[11] V.V. Flambaum, LB. Khriplovich, and O.P. Sushkov, Phys. Lett. B146 (1984) 367. 

[12] V.V. Flambaum, O.P. Sushkov, and V.B. Telitsin, Nucl. Phys. A 444 (1985) 611. 

[13] V.V. Flambaum, LB. Khriplovich and O.P. Sushkov, Nucl. Phys. A 449 (1986) 750. 

[14] LB. Khriplovich, Parity nonconservation in atomic phenomena (Gordon and Breach, Lon- 
don, 1991) 

[15] A.P Platonov, E.E. Saperstein, Nucl. Phys. A486 (1988) 63. 

[16] E.E. Saperstein, M.A. Troitsky, Yadernaja Fizika 1 (1965) 400. 

[17] L.A. Malov, V.G. Soloviev, and I.D. Khristov, Yadernaja Fizika 6 (1967) 1186. 

[18] B. Desplanques, J.F. Donoghue and B.R. Holstein, Ann. Phys. 124 (1979) 449. 

[19] O.P. Sushkov and V.B. Telitsin, Phys. Rev. C 48 (1993) 1069. 

[20] W.C. Haxton, E.M. Henley, and M.J. Musolf, Phys. Rev. Lett. 63 (1989) 949. 

[21] N. Auerbach and B.A. Brown, Phys. Rev. C 60 (1999) 025501. 

[22] V.V. Flambaum, O.K. Vorov, Phys. Rev. C 49 (1994) 1827. 

[23] V.V. Flambaum and G.F. Gribakin, Prog, in Part, and Nucl. Phys., 35 (1995). 

[24] B. Desplanques, Phys. Rep. 297 (1998) 1. 
[25 



S. Tomsovic, Mikkel B. Johnson, A. Hayes, and J.D. Bowman, preprint [nucl-th / 99 1 1 05 7 . 



16 



Table 1: Contributions to the anapole moment of 133 Cs 



K X 10 2 


s.p. 




sv 


Total 










2.5g p + 0Alg n 


5.9g p + 0.4l£ n 


I 


K s X 10 2 


7.0<? p 




1.2g p + 0.24g n 


4.6<? p + 0.24«? n 


II 






-3.5<? p 


1.9g p + 0.17g n 


5Ag p + 0A7g n 


III 






0.% 


-0Ag p - 0.(% n 


-1.2g p - 0.03<? n 


I 


K U X 10 2 


-1.6$, 




-0.3<? p - 0.02<? n 


-1.0g p - 0.02g n 


II 






0.8<? p 


-0.% - 0.02g n 


-1.2g p - 0.02g n 


III 






-0.6# p 


0.9g p + 0.07^ 


-0.2g p + 0.07g n 


I 




-0.5g p 




0.2g p + 0Mg n 


-0.9g p + 0.04<? n 


II 






-0.7g p 


-Q.lg p + 0.02g n 


-1.3g p + 0.02g n 


III 


K c X 10 2 


0.65g pn 






0A5g pn - 0.03g np 


I,H 










0.3Qg pn - 0.02^ 


III 






-3Ag p 


2.9g p + 0A5g n 


4Ag p + 0A5g n + k c 


I 


K tot X 10 2 


4.9g p 




lAg p + 0.25g n 


2.7 g p + 0.25# n + k c 


II 






-3Ag p 


1.5g p + 0A8g n 


2.9g p + 0.18g n + k c 


III 
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Table 2: Contributions to the anapole moment of 205 T/ 



K X 11) 


s.p. 


T/ T/ 


XT / 
V 


total 




^5 


10.8ft, 


-5.7ft, 


4.8ft, + 0.14ft, 


9.9ft, + 0.14ft, 


I 


xlO 2 






3.1ft, + 0.02c/ n 


8.3ft, + 0.02ft, 


II 


Kls 


-2.2ft, 


1-lft; 


-1.0ft, -0.01ft, 


-2.0ft, - 0.01ft, 


I 


X 1U 






— U./ft, — U.Ulft, 


— 1.8ft, — U.Ulft, 


TT 
11 


K-conv 


-0.7ft, 


-0.9ft, 


0.9ft, + 0.23ft, 


-0.8ft, + 0.23ft, 


I 


xlO 2 






-0.6ft, + 0.10ft, 


-2.2ft, + 0.10ft, 


II 


K c X 10 2 


0.85ft, n 






0.64ft, n - 0.06ft,p 


I,H 


K-tot 


7.8ft, 


-5.5ft, 


4.7ft, + 0.35ft, 


7.1ft, + 0.35ft, + k c 


I 


xlO 2 






1.9ft, + 0.1ft, 


4.3ft, + 0.1ft, + k c 


II 



Table 3: Calculated anapole moments for the set of nuclei 



Nuclei 




K x 10 2 


Odd proton nuclei 


133 Cs 


4.9ft, + 0.65ft, n 


2.9ft, + 0.18ft, + 0.36ft,„ - 0.02ft, p 


205 


7.8ft, + 0.85ft, n 


4.3ft, + 0.1ft, + 0.64ft„, - 0.06ft, p 


209 Bi 


5.4ft, + 0.96ft, n 


2.5ft, + 0.3ft, + 0.57ft, n - 0.04ft, p 


Odd neutron nuclei 


135 Ba 


-6.5ft, - 0.25ft, p 


-0.1ft, - 4.6ft, + 0.01ft, n - 0.19ft, p 


137 Ba 


-6.5ft, - 0.25ft, p 


-0.2ft, - 5.7ft, + 0.01ft,„ - 0.23ft, p 


207p b 


-9.6ft, - 0.16ft tp 


-0.1ft, - 6.7ft, + 0.01ft, n - 0.14ft, p 
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